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We explore aspects of the community structures generated by a simple predator-prey model of 
biological coevolution, using large-scale kinetic Monte Carlo simulations. The model accounts for 
interspecies and intraspecies competition for resources, as well as adaptive foraging behavior. It 
produces a metastable low-diversity phase and a stable high-diversity phase. The structures and 
joint indegree-outdegree distributions of the food webs generated in the latter phase are discussed. 
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(S| . I. INTRODUCTION AND MODEL 
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Biological evolution and ecology involve nonlinear interactions between large numbers of units 
and have recently become popular topics among statistical and computational physicists : 1] . How- 
ever, many models used by physicists are unrealistic to the extent of attracting little attention from 
biologists. Here, we introduce a somewhat more realistic model of the dynamics of a predator-prey 
system and explore some aspects of the resulting food- web structures. 
rT"| ■ Recently, we developed simplified models of biological macroevolution [1, Q , in which the repro- 

p I | duction rates in an individual-based population dynamics with nonoverlapping generations provide 

the mechanism for selection between several interacting species. New species enter the community 
through point mutations in a haploid, binary "genome" of length L. The potential species are 
i-Q identified by the index / <E [0, 2 L — 1]. (Typically, only J\f(t) <C 2 L species are present in the 

community at any time t.) At the end of each generation, each individual of species / gives birth 
to a fixed number F of offspring with probability Pj before dying, or dies without offspring with 
probability (1 — Pi). Each offspring may mutate into a different species with a small probability 
fx. Mutation consists in flipping a randomly chosen bit in the genome. 

Here, we consider a model with modified population dynamics that include competition between 
different predators that prey on the same species, as well as a satiation effect for predators with 
abundant prey. Consistent with our previous work Q, the central quantity of the model is an 
antisymmetric interaction matrix M representing predator-prey interactions. Thus, Mtj > and 
Mji < means that / is a predator and J its prey, and vice versa. The elements of the upper 
triangle of M are drawn randomly from a symmetric distribution over [— 1, +1] and kept constant 
during the whole simulation (quenched randomness). A constant, R > 0, represents an external 
resource. The ability of species / to utilize R is rji, which with probability c pro d 1 is chosen to 
be a random number uniform on (0, +1] (i.e., c pro d is the proportion of potential producer species). 
Species with r\i = are consumers. The population size of species I is nr. 

Interspecies competition is modeled by defining the number of individuals of species J that are 
available as prey for J, corrected for competition from other predator species, as 
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where ^P red ( J ) runs over all L such that Mlj > 0. Thus, ^P Icd ( J ) njJ = nj Analogously, 
we define the competition-adjusted external resources available to a producer species I as Ri = 
Rnirji /J2L n L r lL- With these definitions, the total, competition-adjusted resources available to 
species I are 

prcy(/) 

Si = viRi + M u fl u > ( 2 ) 
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where ^P rcy ( / ) runs over all J such that Mjj > 0. 

The functional response of species / with respect to J, is the rate at which an individual of 
species / consumes individuals of J 4, 5]. For ecosystems consisting of a single pair of predator and 
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prey, or a simple chain from a bottom-level producer through intermediate species to a top predator, 
the most common forms of functional response are due to Holling J5|. For more complicated food 
webs, several functional forms have been proposed recently, 0, [a 0, Hi but there is as yet no 
agreement about a standard form. Here, we model intraspecies competition by a ratio-dependent 
[gj Holling Type II 5] form due to Getz [nj, 

Aoj + rij 

where A 6 (0,1] is the metabolic efficiency of converting prey biomass to predator offspring. 
Analogously, the functional response of a producer species / toward the external resource R is 
&m = 'HiRi/i^Si + Tij]. The total consumption rate for an individual of species I is therefore 

^ J Si/m lor AS i < m 

°< = * IR + \ - { l/A forA5z»n I ' 

The birth probability is assumed to be proportional to the consumption rate, Bi = AC/ £ [0, +1], 
while the probability that an individual of / avoids death by predation until attempting to repro- 
duce is 

prcd(7) 

A I = 1- V $ J7 ^. (4) 

j 

The reproduction probability for an individual of species I is Pi(t) = Ai{t)Bj(t). 

As the model is defined above, species forage indiscriminately over all available resources, with 
the output only limited by competition. Also, there is an implication that an individual's total 
foraging effort increases proportionally with the number of species to which it is connected by a 
positive Mjj. A more realistic picture would be that an individual's total foraging effort is constant 
and can either be divided equally, or concentrated on richer resources. This is known as adaptive 
foraging. While one can go to great length devising optimal foraging strategies 0, 0] , we here only 
use a simple scheme, in which individuals of / show a preference for prey species J, based on the 
interactions and population sizes (uncorrected for interspecies competition) and given by 

_ Mijnj 

9lJ — ^ nrov /n ! l°j 

VI R + Ek M IK n K 

and analogously for R by gm — r)iR/[r)iR + J^^k^^ MiKnx]- The total foraging effort is thus 
9iR + J2 P / cy ^ gu = 1< These preference factors are used to modify the reproduction probabilities 
by replacing all occurrences of Mu by Mjjgu and of r\i by r\igin in Eqs. (HJ) — (J3j) . 



II. NUMERICAL RESULTS 



We simulated the model over 2 24 = 16 7 7 7 2 1 6 generations (plus 2 20 generations "warm-up") 
for the following parameters: L = 21 (2 21 = 2 097152 potential species), R = 16 000, F = 2, 
/i = 10~ 3 , Cprod = 0.05, interaction matrix M with connectance C = 0.1 and nonzero elements 
with a symmetric, triangular distribution over [— 1, +1], and A = 1.0. We ran five independent 
runs, each starting from a population of 100 randomly chosen producer species. 

Time series of diversities (effective numbers of species) and population sizes for one run are 
shown in Fig. [IJtop). To filter out noise from low-po pula tion, unsuccessful mutants, the diversity 
is defined as the exponential Shannon- Wiener index [ll|. This is the exponential function of the 
information-theoretical entropy of the population distributions, pi(t) — nj(t)/N to t(t) for the case 
of all species, and analogously for the producers and consumers separately. 

Without adaptive foraging, the system flips randomly between a phase with a diversity near 
ten, and a phase of one or a few producer species with a very low population of many unstable 
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FIG. 1: (Top) Time series of populations (upper curves) and diversities (lower curves) for the adaptive 
model. All species (black), producers (light gray, green online), and consumers (dark gray, red online). 
(Bottom) A representative food web. The producers are shown just above the external resource. Above 
them are the consumers, connected to producers by black arrows and to other consumers by gray (pink 
online) arrows. Arrows point from prey to predator. 



consumer species [12|. Adaptive foraging produces a striking change in the dynamics. There is 
now a metastable low-diversity phase, which gives way at a random time to a stable high-diversity 
phase with much smaller fluctuations. As seen in Fig. QJtop), the switch-over is quite abrupt. 

A representative community food web is shown in Fig. [TJbottom) . This is a "core community," 
extracted from the full community by retaining only specie's with n / > 1 that also existed 25G 
generations earlier. Here, every consumer species preys on at least one producer species, thus 
there are only two trophic levels. (Only links with \Mjj\ > 0.5 are included.) 

Histograms of the correlation coefficient between a species' numbers of prey (indegree) and 
predators (outdegree) are shown in Fig. [^left). The correlations are strongly negative in both 
the simulated full and core communities and also in the majority of the 17 empirical communities 
considered in Ref. H. These negative correlations are explained by the joint indegree-outdegree 
distribution shown in Fig. fright): producers (P) have low average indegree and high average 
outdegree, while consumers (C) show the opposite behavior. 
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FIG. 2: (Left) Correlations between indegree (number of prey) and outdegree (number of predators) 
for full and core communities and 17 empirical communities Q|. (Right) Joint probability distribution 
for indegree and outdegree in core communities. The simulation results in both parts were averaged over 
327 680 communities. 



III. CONCLUSIONS 



We have introduced a model for the biological coevolution of predators and prey, based on the 
ecological concept of functional response. When adaptive foraging is included in the model, it has 
a dynamically stable phase of relatively high diversity. The indegree and outdegree of a species are 
negatively correlated, which is explained by the observation that producers have low indegree and 
high outdegree, while consumers have high indegree and low outdegree. Our model demonstrates 
the high degree of complexity that can be produced, even by simple models of biological evolution 
and ecology. 
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